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PREFACE 


The  United  States  Air  Force  Environmental  Technical  Application  Center 
(USAFETAC)  prepared  this  report  in  answer  to  a  request  from  the  Space  and  Missile 
Test  Center  (SAMTEC) ,  Vandenberg  AFB ,  California,  to  provide  climatological  studies 
of  the  winds  and  temperatures  at  micrometeorological  (micromet)  Towers  300  and  301 
and  to  study  the  feasibility  of  combining  these  data  from  the  two  towers  to  extend 
the  period  of  record  (POR)  of  Tower  301  data.  In  addition  to  a  description  of  the 
programs  utilized  for  the  climatological  summaries,  this  report  gives  an  extensive 
discussion  of  the  regression  procedures  used  to  create  "bogus"  data  for  a  tower  from 
the  data  of  the  other  tower.  It  also  present'*  graphs  of  the  seasonal  comparisons  of 
the  data  from  Towers  300  and  301.  First  Lit  .tenant  Tamzy  J.  Cunningham  and  Major 
Calvin  C.  Naegelin  were  the  original  project  analysts  who  wrote  the  program  designs 
for  the  initial  climatological  summaries. 

If  this  report  is  incorporated  into  another  report  oy  any  agency,  please  give 
USAFETAC  proper  credit  and  furnish  USAFETAC  a  copy  of  the  new  report  if  possible. 
For  further  information  please  contact  USAFETAC. 
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VANDENBERG  AIR  FORCE  BASE 
* WINDS  1  COMPARISON  STUDY 


Wind  Towers  300  and  301 


INTRODUCTION 


The  design  of  the  Space  Transportation  System  (STS)  launch  facility  and  compo¬ 
nents  of  the  space  shuttle  are  sensitive  to  environmental  factors.  The  best  climato¬ 
logical  summaries  of  temperature,  wind  speed,  and  wind  direction  are  essential  for 
the  systems  studies  and  analyses.  The  data  used  in  this  project  are  from  the 
Weather  Information  Network  and  Display  Systems  (WINDS)  at  Vandenberg  AF.O, 
California.  Tower  301,  the  wind  tower  supporting  the  launch  site,  has  a  short 
period  of  record  (POR)  compared  to  that  for  Tower  300.  The  purpose  of  this  project 
was  to  investigate  the  existing  POR,  perform  several  climatological  studies,  and 
determine  the  feasibility  of  expanding  the  POR  for  Tower  301  using  data  from  Tower 
300.  Three  of  the  five  computer  programs  developed  during  this  project  were 
designed  to  answer  specific  problems  for  other  Space  and  Missile  Test  Center 
(SAMTEC)  studies. 


GEOGRAPHY  OF  SITE 

The  STS  launch  facility,  Space  Launch  Complex  6  (SLC-6),  and  micrometeorolog  ical 
(micromet)  WINDS  Tower  301  Elevation  380  ft)  are  located  in  the  southern  foothills  of 
the  Santa  Ynez  Mountains  in  South  Vandenberg's  Punta  De  La  Concepcion  area, 
approximately  1  statute  mile  east  of  Point  Arguello  (Figure  1).  Tranquillon  Ridge 
and  Cypress  Ridge,  with  peaks  ranging  from  1000-1600  feet,  border  the  launch 
facility  and  Tower  301  on  the  east,  northeast,  and  southeast.  Towards  the  north, 
the  south,  and  the  west  the  Santa  Ynez  foothills  slope  into  the  Pacific  Ocean. 
Oceanic  wind  flow  influences  the  tower's  instruments  from  the  east-southeast 
clockwise  to  the  north-northeast.  The  wind  flow  from  the  remaining  quadrant  is 
greatly  influenced  by  the  Santa  Ynez  Mountains. 

WINDS  Tower  300  (elevation  385  ft)  is  located  approximately  4  statute  miles 
north-northeast  of  Tower  301  in  the  northern  foothills  of  the  Santa  Ynez  Mountains. 
The  steep-sided  La  Honda  Canyon,  approximately  800  feet  deep  and  1-1/2  miles  wide, 
is  located  midway  between  Towers  300  and  301.  To  the  north  and  west  of  Tower  300 
the  terrain  gradually  slopes  into  the  ocean;  to  the  northeast  the  terrain  slopes 
into  the  low-lying  area  of  Lompoc  Terrace.  Ridges,  ranging  in  height  from  800-1000 
feet,  border  Tower  300  on  the  south  and  east.  Oceanic  wind  flow  influences  the 
tower  from  the  south-southwestern  to  the  north-northwestern  directions.  Continental 
wind  flow  characterizes  the  remaining  directions. 


DATA  PROBLEMS 

The  micromet  system  is  operated  and  maintained  by  the  Space  and  Missile  Test 
Center  (SAMTEC).  The  data  provide  environmental  information  required  for  planning, 
support  systems  and  facilities.  The  United  States  Air  Force  Environmental  Technical 
Applications  Center  (USAFETAC)  collects,  stores,  retrieves,  and  processes  the 
micromet  data  for  studies  and  analyses  needed  to  support  customer  requirements. 

Tables  1  and  2  show  the  number  of  observations  available  at  Towers  300  and  301, 
respectively.  Data  are  collected  once  per  second  and  averaged  over  a  period  of 
time.  The  basic  time  period  is  30  minutes.  In  the  audit,  an  observation  is 
defined  as  a  date/time  group.  Additional  data  may  be  present  for  one  to  five  levels 
to  include  wind  speed,  wind  direction,  wind  gusts,  and  temperature. 
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Table  2.  Tower  301  Observation  Audit  (Coni' a'* 


Santa  Ynei  Ru«f  Plom 


Figure  1.  Prominent  Terrain  Features,  South  Vandenberg 
{Hinds  and  Nickola,  1968) . 


Tower  300  data  are  available  for  the  period  June  1971  to  May  19 /e.  *lw. 
observations  per  hour  are  available  except  for  a  short  period  from  1  June  197  1  to  / : 
June  1971  when  5-  and  15-minute  observations  were  taken*  May  1978  was  the  la  t  muM!. 
used  in  this  study. 

Tower  301  data  from  February  1966  to  December  1968  were  not  used  in  this  . :  .-jv 
due  to  the  lack  of  temperature  information.  No  data  were  available  froir  1  J  ir.<  >'  ♦ 
until  23  January  1975.  Data  recorded  from  23  January  1975  to  12  April  1975  umt  an 
only  date/time  groups.  After  April  1975  an  observational  count  in  the  60s  indicate' 
the  maximum  data  available.  From  April  1975  through  May  1978  miscellaneous  morrr.s 
are  missing  or  have  low  observation  counts.  May  1978  was  the  termination  point  for 
data  used  in  this  project.  Table  3  gives  a  breakdown  of  the  frequency 
observations  at  the  two  towers. 


Table  3.  Towers  300  and  301  Observation  Frequency. 


TOWER  DATE 


FREQUENCY 


300 

14 

Jun 

71- 

*22 

Jun 

71 

Hour ly , 

5  min 

23 

Jun 

71- 

-31 

May 

78 

Hourly, 

Hour  + 

30 

min 

301 

02 

Feb 

66- 

-31 

May 

67 

Hourly 

01 

Jun 

67- 

-03 

May 

68 

Hourly, 

Hour  + 

10 

min 

04 

May 

68- 

-31 

Dec 

68 

Hourly 

23 

Jan 

75- 

-31 

May 

78 

Hourly, 

Hour  + 

30 

min 

By  studying  the  audit  one  can  see  that  although  there  is  a  lengthy  POR  for  both 
towers,  the  usable  data  are  limited.  The  POR  is  further  decreased  by  some  obser¬ 
vations  containing  only  a  date/time  group  with  no  accompanying  meteorological  data 
or  few  levels.  This  short  POR  is  the  primary  reason  for  attempting  to  supplement 
Tower  301  data  with  data  from  Tower  300. 


PROGRAM  DESCRIPTIONS 

Five  computer  programs  were  developed:  ENBWSVTP,  ENBWTETI ,  ENBWTCMP,  ENBWTCOK , 
and  ENBWTREM .  The  first  two  programs  produce  specific  climatological  data 
summaries.  The  third  program  reformats  the  data  for  use  in  the  fourth  and  fifth 
programs,  which  also  produce  specific  climatoloqical  data  summaries.  The  first  four 
programs  will  be  discussed  in  this  section.  Due  to  the  complexity  and  scope  of  the 
fifth  program,  it  will  be  examined  in  a  separate  section. 

The  cl  ima  tological  studies  were  done  for  five  levels  (12,  54,  102,  204,  and  300 
ft),  the  four  seasons,  and  the  four  wind-directional  sectors  (030°-119°,  120°-209°, 
210°-299°,  300°-029°) .  Each  summary  except  for  the  ENBWSVTP  output  was  computed  for 
daytime  (1400-01002)  and  nighttime  ( 0200-1300Z ) .  The  ENBWSVTP  output  is  summarized 
in  three-hour  blocks.  In  addition  to  the  standard  paper  copies,  microfiche  copies 
are  available.  Due  to  the  volume  of  data  involved,  only  examples  are  shown  in  this 
repor  t . 


Wind  Speed  Versus  Temperature  (ENBWSVTP) 

This  program  assists  in  planning  the  design  of  support  facilities,  systems,  and 
procedures  to  be  used  in  launching  the  space  shuttle.  The  program  calculates 
percentage  frequency  of  occurrence  of  15  wind  speed  and  30  temperature  categories. 
Statistics  were  computed  for  each  of  the  five  levels  for  Towers  300  and  301.  An 
example  of  the  annual  table  for  all  hours  and  wind  directions  for  Tower  301  is  shown 
in  Figure  2.  The  format  of  the  seasonal  table  is  identical.  The  POR  for  Tower  300 
was  14  June  1971-30  September  1977  and  for  Tower  301,  13  February  1975-21  November 
1977. 

External  Tank  Icing  (ENBWTETI) 

The  formation  of  frost  or  ice  on  the  space  shuttle's  external  fuel  tank  is  a 
critical  problem.  The  data  from  this  program  can  be  used  to  determine  the  frequency 
and  severity  of  meteorological  conditions  favorable  to  icing  at  Space  Launch  Complex 
(SLC-6). 
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Figure  2.  Example  of  ENEWSVTP  Output 


The  program  computes  statistical  data  for  wind  speed,  wind  direction,  and  t  - 

ature  using  data  for  a  single  tower.  Two  observations  for  the  same  level,  separator: 
in  time  by  the  user -spec  i  f  ied  interval,  are  located,  and  the  differences  in  wind 
speed,  wind  direction,  and  temperature  for  the  specified  time  interval  are  computed. 
The  difference  is  defined  as  the  later  value  minus  the  earlier  value.  The  mean  and 
the  standard  deviation  of  the  differences  for  all  three  parameters  are  ra 1 cul u* 
tor  the  respective  11  speed  categories.  The  values  are  summarized  for  each  1  *  -  \  *  .  , 
wind  directional  sector,  and  season.  An  al 1 -d irect ional  and  seasonal  summary  is 
also  included.  Figure  3  is  an  example  of  the  ENBWTETI  output.  Only  Tower  301  data 
for  23  January  1975-24  May  1978  were  used. 


Wind  Tower  Comparison  ( ENBWTCMP ) 

The  third  program  provides  an  intermediate  step  to  merge  the  data  from  Towers 
1U0  and  3Ul  together  on  a  tape.  This  output  is  the  input  for  the  comparison  and 
regression  analyses.  The  program  reads  the  POR  tapes  from  Towers  300  and  301  to 
obtain  concurrent  observations  within  the  desired  hour  spans.  If  any  data  ace 
missing,  the  entire  observation  is  discarded.  The  output  is  recorded  on  two 
different  tapes  depending  upon  which  hour  span  the  observation  falls  into.  The 
final  output  tape  contains  concurrent  observations  of  wind  direction,  wind  speed, 
temperature,  and  temperature  differences  for  five  levels  at  each  tower. 


Wind  Tower  Correlation  (ENBWTCOR) 

This  program  compares  the  data  from  Towers  300  and  301.  Using  the  output  from 
the  ENBWTCMP  program,  the  program  compares  wind  direction,  wind  speed,  and  tempera¬ 
ture  from  a  data-rich  tower  (300)  and  a  relatively  data-poor  tower  (301).  Means  and 
standard  deviations  for  differences  in  wind  direction,  wind  speed,  and  temperature 
are  computed  and  categorized  by  11  specified  wind  speeds.  The  difference  is  defined 
as  the  Tower  301  values  minus  the  Tower  300  value.  The  format  of  the  output  table 
is  similar  to  that  of  the  ENBWTETI  data  (Figure  4)  .  The  tables  are  presented 
seasonally  by  wind  directional  sector  and  tower  level.  Also  included  is  an  annual 
and  all-directional  sector  table. 


ENBWTREM 


This  section  describes  the  creation  of  a  new  POR  of  winds  and  temperature  data 
for  micromet  Tower  301.  Due  to  the  limited  amount  of  data  from  Tower  301,  the  ’ata 
from  Tower  300  is  used  to  create  new  "bogus"  data  for  Tower  301. 


Method  of  Analysis 

Initially,  data  from  the  two  towers  were  compared  and  concurrent  observations 
extracted  from  each  wind  tower.  A  statistical  regression  analysis  of  these  compar¬ 
isons  yielded  regression  coefficients  which  were  utilized  to  create  bogus  data  for 
the  data-poor  tower  from  the  data-rich  tower  data. 

The  data  were  divided  into  day  and  night  observations.  The  wind  was  broken  down 
into  U-  and  V-components  for  this  analysis.  The  regression  analysis  compared  the 
five  levels  of  the  U-component,  V-component,  and  temperature  of  Tower  301  against 
the  five  levels  of  Tower  300.  One  component  from  one  level  of  Tower  301  was 
compared  against  the  three  components  of  Tower  300  to  find  the  best  correlation 
(i.e.,  the  U-component  of  level  one  (12  ft)  of  Tower  301  was  compared  against  the 
five  levels  of  the  U-component,  five  levels  of  the  V-component,  and  the  five  levels 
of  the  temperature  component  of  Tower  300)  . 

A  comparison  was  made  between  the  actual  data  of  the  two  towers  before  the  bogus 
data  program  was  implemented.  Besides  being  examined  for  annual  day  and  night 
observations,  the  data  were  also  analyzed  for  seasonal  (spring,  summer,  autumn, 
winter)  categories. 
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Figure  3.  Example  of  ENBWTETI  Output 
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Figure  4.  Example  of  ENBWTCOR  Output. 
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The  statistical  analysis  used  in  this  study  was  a  step-wise  multiple  regression. 
Step-wise  multiple  regression  is  a  statistical  technique  for  analyzing  a  relation¬ 
ship  between  a  dependent  variable  and  a  set  of  independent  variables  (best 
predictors)  in  the  order  of  their  importance.  The  method  used  to  decide  the  ordti 
ot  the  predictors  is  based  upon  the  reduction  of  the  sum  of  squares  of  the  dependent 
variable.  From  this  analysis,  the  regression  coefficients  from  the  best  predictors 
were  used  in  the  creation  of  the  bogus  data.  The  Appendix  to  this  report  provides 
details  of  the  regression  analysis  technique. 

Creating  Bogus  Data 

The  bogus  data  were  created  by  using  the  regression  coefficients  with  the 
formulas 


U(j)  =  UA(J)  +  )  U( 1 }  +  UC(J,i)  V(i ) ]  +  UD(  j )  TD 

i=l 

S 

v(j)  =  VA(j)  +  (VD(J,i)  U(i)  +  VC(J,i)  V( i ) ]  +  VD(j)  TD 
1=1 

*>  5 

T(J)  =  TA(.i)  +  '  [  TP  ( J  ,  1  )  U<  i  )  +  TC(j,l)  V(  1) )  +  X  [F(j,l)  T(i  ) ) 

1=1  1=1 


where  U  =  U-component, 

V  =  V-component, 

T  =  temperature, 

TD  =  temperature  difference  between  level  5  and  level  1, 

j  *  the  level  being  computed,  in  turn,  from  1  through  5, 

i  -  the  level  of  the  second  tower  for  which  the  data  are  available,  also 

for  levels  1  through  5, 

A,  B,  C,  D  and  F  are  the  regression  coefficients. 

The  ENBWTREM  program  uses  the  above  formulas  to  create  a  bogus  POR.  This  pro¬ 
gram  uses  the  regression  coefficients  and  data  values  from  the  data-rich  tower  and 
computes  new  values  for  times  when  no  observations  are  available  for  the  data-  poor 
tower.  If  just  the  data-poor  tower  had  an  observation  or  if  both  towers  had  obser¬ 
vations  at  the  same  time,  the  data-poor  observations  were  used  unchanged.  If  just 

the  data-rich  tower  had  an  observation,  values  were  computed  for  the  data-poor  tower 
using  the  corresponding  date-time  group  from  the  data-rich  tower  observation. 


RESULTS 


Comparison  of  Actual  Data 


The  results  of  the  analysis  are  shown  in  Table  4.  The  table  contains  the  five 
levels  of  each  component.  The  left  column  represents  Tower  301  components  and 
levels.  The  Tower  300  component  that  best  correlates  with  the  corresponding  Tower 
3ol  component  and  level  are  presented  for  day  and  night,  each  season,  and  annual 
summary . 
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ihis  table  shows  that  the  V-component  and  temperature  correlate  better  than  t h«- 
U-eomporient.  Some  of  the  seasons  correlated  better  than  others,  and  the  day 
correlations  were  higheL  than  the  night  correlations.  At  certain  levels  there  is 
not  a  one-to-one  correspondence.  In  some  of  the  U-component  comparisons,  the 
V-component  from  Tower  300  correlated  better  than  the  U-component  at  certain  levels. 
The  results  of  this  analysis  are  also  plotted  (by  component  and  by  time  of  day)  in 
Figures  5-10. 

U-Component  (0200-13002  ,  Figure  5)  .  The  first  and  second  levels  correlate  fairly 
well,  while  there  is~a  large  dif ference  in  levels  3,  4,  and  5.  The  winter,  spring, 
and  fall  pattern  at  levels  3  and  4  do  not  have  very  high  correlations,  while  the 
summer  correlations  are  highest  at  these  levels.  At  level  5,  all  seasons  correlate 
fairly  well,  except  for  summer. 

U-Component  (14QQ-0100Z,  Figure  6).  These  graphs  are  little  changed  from  the  night 
graphs  (Figure  5),  but  exhibit  an  overall  higher  correlation.  The  summer  pattern 
still  has  the  best  correlation,  and  the  winter  the  worst,  except  for  the  first  two 
levels . 

V-Component  (0200-13002,  Figure  7,  and  1400-01002,  Figure  8).  The  pattern  of  the 
V-component  curves  varies  little  between  these  two  f igures .  Overall,  the  summer 
V-component  exhibits  the  highest  correlation. 

Temperature  (020Q-130QZ,  Figure  9,  and  1400-Q1QQZ,  Figure  10).  The  spring  tempera¬ 
ture  curves  show  the  only  significant  variation.  The  best  temperature  correlations 
in  the  spring  are  at  night. 


Comparison  of  Bogus  Data 

Correlation  Analysis.  After  creating  the  bogus  data,  correlation  analyses  were  run 
on  the  bogus  and  actual  tower  data.  The  results  of  these  analyses  are  in  Table  5. 
The  first  two  columns  are  the  correlation  coefficients  for  the  actual  Towers  300  and 
301  data  before  creating  the  bogus  data.  The  next  two  columns  are  the  comparisons 
of  the  actual  Tower  301  data  versus  the  bogus  301  data  created  by  forcing  the  actual 
Tower  301  data  through  the  statistical  analyses. 

There  is  only  a  slight  difference  between  the  correlations  of  the  actual  Towers 
301  and  300  data  and  the  Tower  301  actual  (before)  and  bogus  (after)  data.  The  boqus 
data  have  a  slightly  higher  correlation,  especially  in  the  lower  levels  of  the 
U-component  in  the  night  observations.  The  day  observations  are  also  a  little 
higher.  The  largest  difference  is  in  the  night  U-component. 

Examples  of  Bogus  Data.  Figures  11-14  are  examples  of  the  bogus  data.  The  bogus 
observations  have  no  gust,  integration  period  (IP),  or  reliability  number  (RL)  . 

Figure  11.  The  bogus  data  are  in  the  11G0Z  and  1130Z  observations.  This  is  a 
night  summer  situation  with  light  to  moderate  winds.  The  winds  and  the  surface 
temperatures  correspond,  but  there  is  a  variation  in  the  temperature  differences  at 
upper  levels. 

Figure  12.  This  is  a  night  spring  situation,  and  in  the  bogus  data  there  is 
approximately  a  10-kt  disparity  between  the  actual  wind  and  bogus  winds. 

Figure  13.  This  is  a  day  summer  situation  and  there  is  a  disparity  in  the  wind 
direction  at  all  levels  in  the  17302  bogus  observation.  But  at  2030Z,  the  bogus 
observation  is  not  much  different  from  the  2000Z  actual  observation. 

Figure  14.  In  this  day  winter  situation,  there  is  a  good  correspondence  between 
the  actual  and  bogus  observations  in  both  wind  and  temperature. 


CONCLUSION 


Actual  Data  Comparison 

The  rugged  terrain  at  Towers  300  and  301  is  a  major  factor  in  the  low 
correlation  of  the  U-Component  for  the  towers.  A  contributing  factor  is  the 
prevailing  northerly  wind  component  at  the  two  towers. 
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Figure  5.  Correlation  of  U-Coroponenfc  (02-132)  at  Levels  1-5. 
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Figure  6.  Correlation  of  U-Component  (14-01Z)  at  Levels  1-5 
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Figure  7.  Correlation  of  V-Component  (02-13Z)  at  Levels  1-5. 
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Figure  8,  Correlation  of  V-Component  (14-0 1Z )  at  Levels  1-5 
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Figure  9.  Correlation  of  Temperature  (02-13Z)  at  Levels  1-5. 
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Figure  10.  Correlation  of  Temperature  < 14—012)  at  Levels  1-5 
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Figure  14.  Example  of  Winter  Day  Bogus  Data. 


Table  5-  Comparisons  of  Correlation  Coefficients  for  Towers  300  and  3ul. 


Before  Bogus  After  Bogus 


Component 

Level 

301 (A)/300 (A) 

02-13Z  14-01Z 

301 (A)/301 (B) 
02-13Z  14-01Z 

U 

1 

0.413 

0.634 

0.518 

0.641 

u 

2 

0.445 

0.683 

0.511 

0.687 

u 

3 

0.370 

0.503 

0.560 

0.619 

u 

4 

0.544 

0.651 

0.560 

0.650 

u 

5 

0.562 

0.683 

0.513 

0.689 

V 

1 

0.802 

0.807 

0.860 

0.858 

V 

2 

0.818 

0.813 

0.860 

0.858 

V 

3 

0.821 

0.824 

0.862 

0.858 

V 

4 

0.781 

0.806 

0.829 

0.860 

V 

5 

0.859 

0.862 

0.868 

0.862 

T 

1 

0.859 

0.888 

0.925 

0.885 

T 

2 

0.946 

0.871 

0.947 

0.866 

T 

3 

0.950 

0.873 

0.950 

0.871 

T 

4 

0.849 

0.767 

0.933 

0.871 

T 

5 

0.923 

0.837 

0.925 

0.854 

A  =  Actual  Data 
B  =  Bogus  Data 


CONCLUSION 


Actual  Data  Comparison 

The  rugged  terrain  at  Towers  300  and  301  is  a  major  factor  in  the  lov 

correlation  of  the  U-Component  for  the  towers.  A  contributing  factor  is  the 

prevailing  northerly  wind  component  at  the  two  towers. 

For  each  component  there  is  variance  in  the  level  comparisons.  Again,  the 

terrain  probably  causes  most  of  variances.  The  best  correlations  between  the  two 

towers  occur  in  the  summer  for  all  three  components. 


Bogus  Data  Comparisons 

The  regression  equations  used  did  not  signif icantly  improve  the  correlation 
coefficants  for  bogus  Tower  301 (A) /301 (B)  over  the  coefficants  between  Towers  300 
and  301.  The  examples  in  Table  5  show  that  for  some  levels  the  bogus  values 
correspond  with  the  actual  values.  No  increase  in  accuracy  would  be  achieved  using 
the  bogus  procedure  to  expand  the  POR  of  Tower  301.  The  actual  Tower  300  data  would 
produce  comparable  results. 
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Appendix  A 

REGRESSION  TECHNIQUE 


This  appendix  presents  the  mathematical  equations  that  are  used  by 
ENBWTREM  to  select  the  independent  variables  and  obtain  their  regression  coot i  1  - 
cients.  The  procedure,  which  comes  from  the  IBM  System  360  Scientific  Subroutine 
Package,  Version  III,  starts  with  the  sums  of  cross-products  of  deviations  matrix. 
The  derivation  of  this  matrix  is  not  covered  here;  for  information  on  its  dor  lo¬ 
tion,  refer  to  other  sources  on  regression  analysis. 

In  each  step  of  the  regression 

i  =  1 ,  2,  .  .  .  ,  q 

where  q  =  number  of  independent  variables 

The  first  step  in  selecting  the  independent  variable  entering  the  regression  is 
to  compute  the  amount  of  reduction  of  sum  of  squares  for  each  variable  by 


where 


a  .  . 
ID 


Cj  =  reduction  of  sum  of  squares, 

ajj  =  an  element  in  the  sum  of  cross-products  of  deviations 
matrix  which  will  be  modified  in  successive  steps, 

j  =  1,  2,  .  ..,  q  are  independent  variables  (j  /  variables  deleted 
and  variables  entered  before  the  ith  step), 

y  *  dependent  variable. 

The  second  step  is  to  find  the  maximum  (over  j)  of  C. 

Set  Si  =  Cj  to  indicate  the  sum  of  squares  that  will  be  reduced  in  the 

ith  step.  The  proportion  of  Sy  to  the  total  is  obtained  by 


D 


n 

v  _  p 

where  D  *  _  (yj-y)  ' 

J-l 

n  =  number  of  observations. 

The  cumulative  sum  of  squares  reduced  is  obtained  by 

Scum  =  Scum  +  Si  [&-?'> 

and  the  cumulative  proportion  reduced  by 

pcum  =  pcum  +  P  (A-M 

The  multiple  correlation  coefficient  is  computed  by 

R  -  (Pcum) 1/2 


and  adjusted  for  degrees  of  freedom  by 


The  F-value  for  analysis  of  variance  is  given  by 

^  k 
“  cun. 

(l)-Scum)  /  (n-k-1) 

The  standard  error  of  the  estimated  y  is  obtained  by  the  use  of  the  formula 


y .  1  . .  i 

and  adjusted  by 

sc  -  s  ( (n-i ) / (n-k ) j 
Then,  the  following  is  computed 


,0-:'  vl/? 

«  ! _ -U-M 

n-k-f  J 


1/2 


a  -  ■  *  a.  .  + 


til 


33  33 


(A- 


where  i  =  variable  entered  in  the  i tb  step, 


where 


3  =  v, ,  v 


1'  2'  *  V /i-1 


before  the  itn  step,  and 


are  the  variables  entered  in  the  regression 


a,k 

a  . ,  -  —  • 

-ik  au 


k  *  1,  2,  .  ..,  m  are  variables  including  y  (k  /  variables  deleted  and  the 
variables  entered  in  the  i step). 


Regression  coefficients  are  computed  by: 


b  .  ss  g  . 

i  Jiy 


k)i-l  ^(i-l)y  ”  biq(i-l)i 


b. 


2  =  q(i-2)v  ~  biq(i-2li  ~  bi- lq ( i-2) ( i-1 ) 


etc.  , 

and  the  value  of  the  intercept  is 


f  /, .  1  ;>  'j 


bo  =  V  '  l  b.i. 

]*1 


(A -Jl) 


where  k  =  number  of  independent  variables  in  the  regression. 

Standard  error  of  regression  coefficients  are  given  by 

sb j  =  (aj jsy.l2. . . i* X/^  (A-l£) 

where  j  =  v^,  v2,  ...,  v^  are  variables  i,  the  regression,  and  t-values  as 


J  s 


bj 

Perform  the  reduction  to  eliminate  the  variable  entered  in  the  i^h  step, 

a  ,  *  a  _  a  a , 

3k  jk  ji  ik 


(A-131 


(a-im 
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